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Abstract 

We present an exact solution of the equations for orbit determi- 
nation of a two body system in a hyperbolic or parabolic motion. 
In solving this problem, we extend the method employed by Asada, 
Akasaka and Kasai (AAK) for a binary system in an elliptic orbit. 
The solutions applicable to each of elliptic, hyperbolic and parabolic 
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orbits are obtained by the new approach, and they are all expressed in 
an explicit form, remarkably, only in terms of elementary functions. 
We show also that the solutions for an open orbit are recovered by 
making a suitable transformation of the AAK solution for an elliptic 
case. 

keywords: astrometry — celestial mechanics — orbit determination 

1 Introduction 

Two body problems are very classical in celestial mechanics and have been 
studied thoroughly since Kepler discovered the laws of motion of celestial 
objects (e.g., Aitken 1964, Goldstein 1980, Danby 1988, Roy 1988, Murray 
and Dermott 1999, Beutler 2004). The regular orbits of a system of two 
masses in Newtonian mechanics are of three types: ellipse, parabola and 
hyperbola. The latter two cases, in which the separation of two bodies will 
become infinite in the remote future, may be called open orbits. The singular 
orbit is a linear one corresponding to a head-on collision, which is extremely 
special. In this paper, we consider regular orbits mentioned above: for a 
binary system, the orbit determination may bring us some informations about 
its formation and evolution mechanism. In an open orbit, one may infer an 
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impact parameter and an initial relative velocity of two masses, one of which 
may be ejected by some explosive mechanism such as a supernova or by three 
body scattering. For instance, some observations reveal that kicked pulsars 
move at unusually high speed (Anderson et al. 1975, Hobbs et al. 2005). 

The orbit determination of visual double stars was solved first by Savary 
in 1827, secondly by Encke 1832, thirdly by Herschel 1833 and by many 
authors including Kowalsky, Thiele and Innes (Aitken 1964 for a review on 
earher works; for the state-of-the-art techniques, e.g, Eichhorn and Xu 1990, 
Catovic and Olevic 1992, Olevic and Cvetkovic 2004). Here, a visual binary 
is a system of two stars both of which can be seen. The relative vector from 
the primary star to the secondary is in an elliptic motion with a focus at 
the primary. This relative vector is observable because the two stars are 
seen. On the other hand, an astrometric binary is a system of two objects 
where one object can be seen but the other cannot like a black hole or a 
very dim star. In this case, it is impossible to directly measure the relative 
vector connecting the two objects, because one end of the separation of the 
binary, namely the secondary, cannot be seen. The measures are made in 
the position of the primary with respect to unrelated reference objects (e.g., 
a quasar) whose proper motion is either negligible or known. 
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As a method to determine the orbital elements of a binary, an analytic 
solution in an explicit form has been found by Asada, Akasaka and Kasai 
(2004, henceforth AAK). This solution is given in a closed form by requiring 
neither iterative nor numerical methods. One may naturally seek an analytic 
method of orbit determination for open orbits. An extension for open orbits 
done earher by Dommanget (1978) used the Thiele-Innes method, namely 
solved numerically the Kepler equation. As a result, the method by Dom- 
manget does not provide an explicit solution in a closed form. Therefore, let 
us extend the explicit solution by AAK to open orbits. Then, we would face 
the following problem. AAK formalism uses a fact that the semimajor and 
semiminor axes of an ellipse divide it into quarters. This fact plays a crucial 
role in determining the position of the common center of mass on a celestial 
sphere; we should note here that the projected common center of mass is 
not necessarily a focus of an apparent ellipse. The division into quarters is 
possible for neither a parabola nor a hyperbola. 

The purpose of this paper is to generalize AAK approach so that we can 
treat an open orbit. This paper is organized as follows. Sec. 2 presents a 
generalized AAK formalism. In Sec. 3, the generalized approach is employed 
to obtain the method of orbit determination for a hyperbolic orbit. In Sec. 
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4, the formula for a parabolic orbit is presented. In Sec. 5, we recover these 
formulae by making a suitable transformation of that for an elliptic orbit 
with some limiting procedures. Sec. 6 is devoted to Conclusion. 



2 Generalizing AAK formalism for an elliptic 
motion 

2.1 An apparent ellipse 

We denote by {x, y) the Cartesian coordinates on a celestial sphere that is 
perpendicular to the line of sight. A general form of an ellipse on a celestial 
sphere is 



which is characterized by five parameters; the position of its center, the length 
of its semimajor/semiminor axes and the rotational degree of freedom. By 
at least five measurements of the location of a star, one can determine all 
the parameters. Henceforth, we adopt the Cartesian coordinates (x, y) such 
that the apparent ellipse can be reexpressed in the standard form as 



ax"^ + 13 f + 2'-^xy + 25x + 2sy=l, 



(1) 




(2) 
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where we assume a > b without loss of generahty. The elhpticity, e, is 



2.2 The common center of mass 

A focus of the original Keplerian ellipse is not always that of the apparent one 
because of the inclination of the orbital plane. We should note that a focus of 
the original Keplerian orbit is the common center of mass of a binary, around 
which a component star moves at the constant-areal velocity following the 
Keplerian second law (the conservation law of the angular momentum in the 
classical mechanics). This enables us to find out the location of the common 
center of mass as shown below. 

A star is located at Pj = {xj,yj) = {a cos u j,b sin uj) on a celestial sphere 
at each epoch tj for j = 1, ■ ■ ■ , 4, where tj > tk for j > k. Here, Uj denotes 
the eccentric angle in the apparent ellipse but not the eccentric anomaly in 
the true one; the eccentric anomaly of the original Keplerian orbit is not 
observable. We assume anti-clockwise motion, such that uj > Uk for j > k. 
All we must do in the case of the clockwise motion is to change the signature 
of the area in Eq. (0)) in the following. We define the time interval as 




^) = - tk- 
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The common center of mass of the binary is projected onto the celestial 
sphere at Pq = {xc, yc)- Even after the projection, the law of constant-areal 
velocity still holds, where we should note that the area is swept by the line 
interval between the projected common center and the star. The area swept 
during the time interval, T{j,k), is denoted by S{j,k). The total area of 
the observed ellipse is denoted by S* = nab. The law of the constant areal 
velocity on the celestial sphere becomes 

S Sij,k) 



T T{j, k) ' 

where 



(3) 



Sij, k) = -ah Uj — Uk -ismuj — sinuk) + ^(cosm,- — coswfc) . (4) 

2 ^ a '1 ^ ' 

Equation (jH)) is rewritten explicitly as 

^(2,1) S(3,2) 



r(2,l) T(3,2)' 

SM ^ 5(4,3) 
T(3,2) T(4,3)' 

They are solved for xq and yc as 



(5) 
(6) 



B1C2 — B2C1 , , 



AIB2 - A2BI 

^ CrA2-C2A, 



where 

= T{j + 1, j) sin Uj+2 + T{j + 2, j + 1) sin uj 

-T{j + 2,j)smuj+i, (9) 

Bj = T{j + 1, j) cos Uj+2 + T{j + 2, j + 1) cos Uj 

-TU + 2J)cosUj+i, (10) 

Cj = T{j + l,j)u,+2 + Tij + 2,j + l)uj-T{j + 2,j)u,+i. (11) 

The periastron is projected onto the observed elhpse at = (x^, ypj = 
(a cosM^, bsinujij. The ratio of the semimajor axis to the distance between 
the center and the focus of the elhpse remains unchanged, even after the 
projection. Hence, we find 

Pa = ;^Pc- (12) 

The positional vector is still located on the apparent ellipse given by Eq. 
Q. We thus obtain the ellipticity as 



^K = l/# + i- (13) 



2.3 Projection onto a celestial sphere 

In the original derivation of AAK formula, the fact that the semimajor and 
semiminor axes divide the area of the ellipse in quarters. This still holds 



even after the projection. Namely the projected semimajor and semiminor 
axes divide the area of the apparent elhpse in quarters, though the original 
semimajor and semiminor axes are not always projected onto the apparent 
semimajor and semiminor ones. This way of the derivation, however, can be 
used for neither hyperbolic nor parabolic cases, where there are no counter- 
parts of the semiminor axis. Hence we shall employ another method. 

In this paragraph, we use the Cartesian coordinates (X, Y) on the original 
orbital plane. Let i be the inclination angle between the original orbital 
plane and the celestial sphere. We define as u; the angular distance of the 
periastron, namely the angle between the periastron and the ascending node. 
Let us express the original Keplerian ellipse as 

^ + ^ = 1- (14) 
We consider the line that is perpendicular to the semimajor axis at a focus 
(^K^K' '-')■ This line intersects the original ellipse at points Q — (^k^K' ^k(1" 
ej^)) and R = (^k^K' ~'^k(-'- ~ ^k))' later convenience, we adopt the 
coordinates {X,Y) whose origin is located at the focus, by making a trans- 
lation as X = X — ox^K- Then, one rewrites Q = (0, aj^(l — Cj^)) and 
R=(0,-aK(l-e|.)). 

Only in this paragraph, we adopt other Cartesian coordinates {x',y') so 



that the ascending node can be located on the x'-axis and the origin can be 
the common center of mass. The true periastron of the original ellipse is pro- 
jected at = Wp^)y'p^ = a}<^{l — eYj{cosu!,smu!COsi). The point denoted 
by Q is projected at Pq = (xqj^q) = (^K^^ ~ ej^)(— sina;,cosa;cosi). It 
is useful to consider the following invariants because the components of a 
vector depend on the adopted coordinates. 

|P^|2 = aj^{l - ej^)^(cos^a; + sin^a;cos^i), (15) 
|Pq|2 = aj^{l - e|^)^(sin^u; + cos^ a; cos^ i) , (16) 
|Pa X PqI = «k(l + eK)(l - ^k)' cosz. (17) 
We consider the area surrounded by the ellipse and the line interval be- 
tween Q and R. This area is divided into equal halves by the semimajor axis. 
Even after the projection, the divided areas are still equal halves. Hence one 
can determine the location of the projected Q as 

in the apparent ellipse coordinates, where we defined 

s^^ ^. (19) 

In this computation, it is useful to stretch the apparent ellipse along its 
semiminor axis by a/b so that one can consider a circle with radius a. In 
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this stretching, importantly, the areal division into equal halves still holds. 



We make a translation as x ^ x — xc and y y — yc so that the projected 
common center of mass can become the origin of the new coordinates. Then, 



we have 



Pa 



1 - ex 



aycs hxcs 



(20) 
(21) 



Hence, we obtain the invariants from these vectors as 



1 -e 



K 



(4 + l/c)> 



{aW^ + b'xl){l-el) 

Pr' 



ipa X PqI = «Ki-eK)v'i-e; 



^2 

K' 



(22) 

(23) 
(24) 



whose values can be estimated because a, 6, xc, yc and ej^ have been already 
all determined up to this point. 

Equations (fT5|l - (fT7jl for cosz, aj^ and cos 2^7 are solved as 



cos I 
cos 2uj 



1 



e-4), 



1 - e 



K 



\ 



:i+eK)'|PAP+|PQ 

1 + cos^ i 



:i + *k)'|PaP-IPqI 



e|^)^ sin^ i 



(25) 
(26) 
(27) 
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where we define 

^ (l + eK)W + |PQP 
^ (l + eK)|PAxQQl ■ ^ ' 

One can show 

e > 2, (29) 

because the arithmetic mean is not smaller than the geometric one. It is 
worthwhile to mention that Eq. (|25|) is obtained by solving a quadratic 
equation for cosi as 

cos^ i — i cos i + 1 = 0, (30) 

which can be obtained from Eqs. ()15j) - (fT7j) by eliminating uj and aj^. Fur- 
thermore, one can prove that a root of cosi = (^ + ^J^^ — A)/2 must be 
abandoned because Eq. (j29p implies that it is always larger than the unity. 
Only in the case of z = 0, the apparent ellipse coincides with the true orbit. 
Hence, the ascending node and consequently the angular distance of the pe- 
riastron make no sense. As a result, the denominator of R. H. S. of Eq. (j27|) 
vanishes. 

Equations (UHl), (ESI), and (j27|) agree with those of AAK, where dif- 
ferent notations were employed. In this paper, the semiminor axis is not used 
for areal divisions. Therefore, this formalism can be generalized straightfor- 
wardly to an open orbit, as shown below. 
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3 The solution for a hyperbolic orbit 



3.1 An apparent hyperbola 

Let a star move in a hyperbola on a celestial sphere. Without loss of gener- 
ality, we can assume that the hyperbola is expressed as 

^-^ = 1' (31) 

and the orbit is the left-hand side of the hyperbola, x < 0. Then the position 

of the star at each epoch tj is denoted by 

Pj ~ {■^j-iVj) 

— (— a coshtij, 6sinhMj). (32) 

The projected common center of mass Pq = (xq,^^) is not necessarily 
a focus of the apparent hyperbola but the projected focus of the original 
Keplerian hyperbola. 

The projected areal velocity with respect to the projected common center 
of mass is denoted by dS/dt. The law of the constant areal velocity on the 
observed plane is written as 

dl _ S{j, k) 

dt-mky ^''^ 
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where for tj > tk we obtain 



1 



'SU, k) = —-ah Uj — Uk -\ (sinhuj — sinhuk) + — (coshwj — coshwfc) .(34) 

^ a 

Equation (jH^ is rewritten explicitly as 

^(2,1) 3(3,2) 



r(2,i) r(3,2)' 

g(3,2) ^ g(4,3) 

T(3,2) T(4,3)' 

They are solved for xq and yc as 



D1E2 — D2E1 ' 
'Z^i^s - D2E1 

where 



= T(j + 1, j) cosh + r(j + 2, j + 1) cosh 



u 



(35) 
(36) 



E1F2 — E2F1 . 
= a-— -— , (37) 



FJD2-F2D, 
yc = b-— — — , (38) 



Dj = T{j + l,j)smh.Uj+2 + T{j + 2,j + l)smh.Uj 

-T{j + 2,j)smhuj+i, (39) 



-T(j + 2,j)coshM,+i, (40) 
F, = T{j + l,j)u,+2 + Tij + 2,j + l)u,-Tij + 2,j)u,+i. (41) 

The periastron is projected onto the observed hyperbola at = (x^, ypj. 
The ratio of the semimajor axis to the distance between the center and the 
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focus of the hyperbola remains the same, even after the projection. Hence, 
we find 

Pa = ^Pc- (42) 

The positional vector is still located on the apparent hyperbola given by 
Eq. (jnH). We thus obtain the ellipticity as 



- - i- («) 



3.2 Projection onto a celestial sphere 

In this paragraph, we use the Cartesian coordinates (X, Y) on the original 
orbital plane. Let us express an original Keplerian hyperbola as 

^-^ = 1- (44) 

We consider the line that is perpendicular to the semimajor axis at a focus 
('^K^K' This line intersects the original hyperbola at Q = (ax^K' ^K^^K~ 
1)) and R = (ctK^K' ~^K(^k ~ 

Only in this paragraph, we shall employ other Cartesian coordinates 
{x', y') so that the ascending node can be located on the x'-axis and the 
origin can be the common center of mass. The true periastron of the original 
hyperbola is projected at = {^'p^iV'p^ = '^k(^K l)(cosu;, sinu; cosi), 
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where i and cu are the inchnation angle and the angular distance of the perias- 
tron, respectively. The point denoted by Q is projected at Pq = (^^^q, Uq) = 
'^k(^K ~ sino;, cosa;cosi). We shall use the following invariants as 

I-^aI'^ = Oq^{e-^ — l)^(cos^ u + sin^ u) cos^ i), (45) 
|Pq|2 = aj^(ej^ - l)^(sin^u; + cos^u;cos^i), (46) 
|Pa X PqI = «k(eK + l)(eK - l)'cosi. (47) 

We consider the area surrounded by the hyperbola and the line interval 
between Q and R. This area is divided into equal halves by the semimajor 
axis. Even after the projection, the divided areas are still equal halves. Hence 
one can determine the location of the projected Q as 

PQ-^^c--^,yc + ^J, (48) 
in the apparent hyperbola coordinates, where we defined 



^h - ^ ■ (49) 

We make a translation as x ^ x — xc and y ^ y — yc so that the center of 
the coordinates can move to the projected common center of mass. Then, 
we have 

PA^^—^i^cyc), (50) 
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where we used Eqs. (jl^ and PHj) . Hence, we obtain the invariants from 
these vectors as 



Al 



- 1 



a'yl + b'xl){e'^-l) 



Qi 



|PAxPQ| = a6(eK- 1)^4-1 



(52) 

(53) 
(54) 



whose values can be estimated because a, b, xc, yc and e-^ have been all 
determined up to this point. 



Equations (j45|l - (j47|l for cosi, aj^ and cos 2a; are solved as 



cosz 
cos 2uj 



2(^-ve-4), 



K 



(gK + l)^|PAl^ + |PQl 



(eK + l)'|PAl 



1 + cos^ i 



where we define 



(eK + l)'|PAP + |PQl 



(eK + l)|PAxQQl 
In the similar manner to the elliptic case, one can show 

^ > 2. 
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(55) 
(56) 
(57) 



(58) 



(59) 



Hence, for a quadratic equation for cosi as 

cos^ i — ^ cos i + 1 = 0, (60) 

which is derived from Eqs. ()45|) - (jT7j) . One can prove that a root of cosi = 
+ VC^ ~ 4) /2 is larger than the unity according to Eq. (IHTIjl and thus must 
be abandoned. 



4 The solution for a parabolic orbit 

4.1 An apparent parabola 

Let a star move in a parabola on a celestial sphere as 

+ 4qx = 0. (61) 

Then the position of the star at each epoch tj is denoted by 

Pj = (xj,yj) 

= {-^qulV2qu,). (62) 

The projected common center of mass Pq = {xQ,yQ) is not necessarily 
a focus of the apparent parabola but the projected focus of the original 
Keplerian parabola. 
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The law of the constant areal velocity on the observed plane is written as 

dS S{j, k) 



dt T{j,k)- 



(63) 



where for tj > tk we obtain 



SU, k) 



+ - 



3 
3 



-X 



j + V-Xk)yc + ^Vq^c 



Equation ()(i3|l is rewritten explicitly as 

5(2,1) 3(3,2) 



r(2,l) 
T(3,2) 

They are solved for xc and yc as 



r(3,2)' 

^(4,3) 

T(4,3)- 



yc 



Hih — H2I1 
G1H2 — G2H1 ' 
G1I2 — G2I1 



G1H2 — G2H1 



where 



(64) 



(65) 
(66) 



(67) 
(68) 



3[T(j + 1, jV-x,+2 + T(j + 2, J + 1) 



-T{j + 2,j)^-Xj+i 
3 



(69) 



[T{j + 1, j>,+2 + T(j + 2, J + l)x, - T{j + 2, j>,+i], (70) 



-[T{j + l,j)xj+2J-Xj+2 + T{j + 2,j + l)xj 



-TU + 2,j)xj+iJ-Xj+i 



(71) 
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The periastron is projected onto the observed parabola at = {xp^, ypj . 
The semimajor axis is projected onto a hne, which may be expressed as 
y — Kx + L. This hne intersects the apparent parabola only at the pro- 
jected periastron. Therefore, we find K — Q because \x\ must be larger than 
^J—x for a sufficient large \x\. In addition, the projected semimajor axis 
goes through the projected common center. This implies L — y^. Hence we 
obtain 

Pa-(-|.c). (72) 
4.2 Projection onto a celestial sphere 

In this paragraph, we use the Cartesian coordinates {X, Y) on the original 
orbital plane. Let a Keplerian parabola be 

y2 + Aqy^X = 0. (73) 

We consider the line that is perpendicular to the semimajor axis at a focus 
(~9K' '^^^^ intersects the original parabola at Q = (— ^^k) 

R= (-?K'-2?k)- 

Only in this paragraph, we employ other Cartesian coordinates {x',y') 
so that the ascending node can be located on the x'-axis and the origin 
can be the common center of mass. The true periastron of the original 
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parabola is projected at = (^^)?/^) — Q'j^(cosa;,sina;cosi), where i 
and LU are the inchnation angle and the angular distance of the periastron, 
respectively. The point denoted by Q is projected at Pq = {x'q^i/'q) = 
2gj^(— sino;, cosa;cosi). We shall use the following invariants as 

|P^|^ = g|^(cos^a; + sin^a;cos^i), (74) 
|Pq|^ = 4g|^(sin^a; + cos^cucos^i), (75) 

|PaxPqI = 25kcos^- (76) 
We consider the area surrounded by the parabola and the line interval 
between Q and R. This area is divided into equal halves by the semimajor 
axis. Even after the projection, the divided areas are still equal halves. Hence 
one can determine the location of the projected Q as 

Pq= (^c-^,yc + «p), (77) 

in the apparent parabola coordinates, where we defined 

sp = V-(^c + 4gxc). (78) 

We make a translation as x x — xc and y y — yc so that the origin 
of the coordinates can be the projected common center of mass. Then, we 
obtain 
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ycfp 



(80) 



Hence, we obtain the invariants from these vectors as 

,2 {yl + 4qxcY 



A I 



IPaxPqI 



16g2 

(yg + 4gxc)(yg + 4g^) 

4g2 

[-(yi + 4ga;c)]3/^ 



4g 



(81) 
(82) 
(83) 



whose values can be estimated because q, Xq and Uq have been all determined 
up to this point. 

Equations (f7I|l - (ffn|l for cosi, a-^ and cos2cij are solved as 



cos I 

IK 
cos 2uj 



(^-\/e-4), 



4|PaP + |PqI 

2 ^ 1 + cos2 i 

4|PaP-|PqP 

4gj^ sm I 



where we define 



4|PaI^ + |PqP 
2|PaxQq| • 

In the similar manner to the above two cases, one can show 



(84) 
(85) 
(86) 

(87) 



e>2. 



(88) 



Hence, for a quadratic equation for cos i as 



cos i — cos 2 + 1 = 0, 



(89) 
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which is derived from Eqs. (j74|) -(|76| ) . One can prove that a root of cosi = 
+ V'C^ ~ 4)/2 is larger than the unity according to Eq. (jHH|) and thus must 
be abandoned. 



5 Transformations from an elliptic case to hy- 
perbolic/parabolic ones 

5.1 To a hyperbolic case 

Let us rederive the formula for a hyperbolic case from that for an elliptic one 
by making a transformation as 

u iu, (90) 

a — > —a, (91) 

b -lb, (92) 

which imply 

cosM — >• coshii, (93) 

sin M i sinh u, (94) 
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where i = y/^. Then, from Eqs. Q-dlll) and ((23)- (SD) we find 

A, ^ Id,, (95) 

B, ^ E„ (96) 

C, ^ iF,. (97) 

We can thus show that the location of the common center is transformed 
from Eqs. ((Tj) and (jH)) to Eqs. (jH7|) and (jHHjl . Equation (fT5|) is transformed 
into Eq. (gS)), Eqs. into Eqs. (jS21)-(Ell), Eqs. (ESI)- (EH) into Eqs. 

(j^ -(|37 j) . because ^ remains unchanged. 

5.2 To a parabolic case 

To rederive the formula for a parabolic case, we perform a transformation 
from an elliptic case with a limiting procedure as 

e ^ 1, (98) 
q = lima(l-e), (99) 
x' = x-a, (100) 

where the finite q implies a — > oo and 6^ = a(l + e) x a(l — e) ^ 2ag. Then, 
we find 

A- - (101) 
24 



B, Ih,, (102) 

- + (103) 

We can transform the location of the common center from Eqs. ((7)) and ((HI) 

to 

H1I2-H2I1 , - 

l/c ^ -V^TT-f} (105) 

which agrees with Eq. (jUHjl . We thus recover Eq. as 

4 - (106) 
Equation (fT^ is transformed as 

^ 1, (107) 



where we used a 00 and 6^ — 2aq. By using Eq. ()1()7|1 and 6^ — 2ag, we 
obtain 

IPaP ^ (108) 



IP ,2 ^ 2[a^yg + 4aV(a + x-e)^](l-eK) 

(yg + V)(yg + 4gx^c) 
^ 4g2 
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(109) 



IPaxPqI 



4q 



;iio) 



which agree with Eqs. (jRT]) - !^^ . 

Equations (j^^ -()27 |) are transformed as 



cos? 



cos 2u 



2(e-v^'-4), 



hm ai^fl — e-fT-] 



hm \ 



:i + eK)^|PAP+|PQP 
1 + cos^ i 



-2\ 



*|PaI' 



1 + cos^ i 

:i+«KnPAp 



IPqI'^ 



qI 



4|PaP 



lPr'2 



4g|r- sin^ i 



;iii) 



;ii2) 



;ii3) 



where ^ remains unchanged. They agree with Eqs. (j84p - (j86p . 



6 Conclusion 

The formulae for orbit determination of ehiptic, hyperbohc and parabohc 
orbits are obtained in a unified manner by generahzing AAK approach, which 
originally needed a fact of the areal divisions by the semimajor and semiminor 
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axes of an ellipse. We show also that the present formulae are recovered from 
AAK result by a suitable transformation among an ellipse, hyperbola and 
parabola. 
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